Delta-Function Potential with a Complex Coupling 



Ali Mostafazadeh 



Department of Mathematics, Koq University, 
34450 Sariyer, Istanbul, Turkey 
amostafazadeh@ku . edu.tr 



Abstract 

We explore the Hamiltonian operator H = —-^ + z 5{x) where x G M, 5{x) is the Dirac 
delta function, and z is an arbitrary complex coupling constant. For a purely imaginary 
z, H has a spectral singularity ai E = —z^ /A S For 3^(2;) < 0, H has an eigenvalue 
ai E = — For the case that ^{z) > 0, H has a real, positive, continuous spectrum 
that is free from spectral singularities. For this latter case, we construct an associated 
biorthonormal system and use it to perform a perturbative calculation of a positive-definite 
inner product that renders H self-adjoint. This allows us to address the intriguing question of 
the nonlocal aspects of the equivalent Hermitian Hamiltonian for the system. In particular, 
we compute the energy expectation values for various Gaussian wave packets to show that 
the non-Hermiticity effect diminishes rapidly outside an effective interaction region. 

PACS number: 03.65.-w 

Keywords: Complex potential, delta function, pseudo-Hermitian, inner product, metric op- 
erator, spectral singularity, PT-symmetry 
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1 Introduction 



The observation that a complex potential can define a consistent unitary quantum system has 
recently led to a considerable research activity. This is mostly focused on the study of the complex- 
valued PT-symmetric potentials v for which the Hamiltonian 



has a real discrete spectrum. In general the eigenvalue problem for H is defined along an ap- 
propriate contour F in the complex plane with suitable boundary conditions at infinity pj. This 
allows one to identify H with a densely defined and generally non-self-adjoint operator acting in 
a separable Hilbert space ?i, with the typical choice for H. being L'^{T). 

The problem of whether and how one can formulate a consistent quantum system having 
H as its Hamiltonian has found a satisfactory solution within the context of pseudo-Hermitian 
quantum mechanics [21 [3l HI [5] . It turns out that H must be diagonalizable. In particular, there 
must exist a complete basis of Ti consisting of eigenvectors of H . This is a physical requirement 
of the standard quantum measurement postulate [21 E]- For a diagonalizable Hamiltonian with a 
discrete spectrum, the reality of the spectrum is equivalent to the existence of a positive-definite 
inner product (■,■)+ that renders H self-adjoint [3, |8]. The latter is a necessary and sufficient 
condition for the existence of an equivalent Hermitian Hamiltonian h that acts in H, [3 [9]. These 
observations lead to the realization that the physical system under investigation may be equally 
well described by a Hermitian Hamiltonian within the framework of the conventional quantum 
mechanics [IOl[IIl[2]. 

The above discussion also applies to complex potentials that are not "PT-symmetric. The 
purpose of this paper is to study in detail one of the simplest (though highly nontrivial) examples 
of complex potentials that happens not to be PT-symmetric, namely the delta-function potential 
with a complex coupling: 



We wish to explore the possibility of defining a unitary quantum system based on the standard 
Hamiltonian ([T|) and the potential ([2|). The reference Hilbert space [2] is given by ?^ = L^(M) 
and it is not difficult to see that the spectrum of the Hamiltonian fails to be discrete. Because 
of the presence of the continuous part of the spectrum the results reported in [121 Ej may not 
hold. Nevertheless, they provide some useful guiding principles that we will follow. Specifically, 
we will attempt to construct a biorthonormal system whenever possible and use it to define an 
appropriate positive-definite inner product that renders H self-adjoint. An example of a successful 
application of this strategy is the scattering potential [S] 



H 




(1) 



C G C. 



(2) 




for |a;| < ^ 
for \x\ > "I, 



(3) 
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where A G M, L G M+, and 

{—1 for a; < 
for a; = 
1 for X > 0. 

Note that unhke (l3l), the delta- function potential ([2]) fails to be "PT-symmetric, and as we will 
see, depending on the value of (, it may lead to the presence of a spectral singularity |I3j. 

Before starting our analysis of the properties of (E]), we wish to point out that complex poten- 
tials consisting of one or more delta functions have been studied in [H] - [21] and that the issue of 
the emergence of spectral singularities for PT-symmetric potentials has been considered in |22j . 



2 Spectral Properties and Biorthonormal Systems 

Consider the time-independent Schrodinger equation, 

Hijj = Ei/j (4) 

subject to bounded boundary conditions at x = ±oo. Substituting ([T]) and (E]) in (jll) and intro- 
ducing the dimensionless quantities: 

X 2miC ^ 2me'^E 

where £ is an arbitrary length scale, we can express (jl]) in the form 

- z^"(x) + z5{x)iIj{x) = Ez^(x). (6) 
Clearly E belongs to the spectrum of the dimensionless Hamiltonian 

H:=^i^ = -^ + ^^(x). (7) 
The solution of ([6]) has the form 

V,(X)=*(X):=P-C^''-C W 

^ ^ ^'^^ ^ I A+e*^" + E+e-*'^" for x>0, ^ ^ 

where k := V^, the coefficients A^,B^ G C are arbitrary but not both vanishing, i.e., + 
^ 0, and 

= (1 - ^)A_ - ^ fi+ = ^ A_ + (1 + (9) 

^ ^ 2k' 2k ' ^ 2k ^ 2k' ^ ' 

A straightforward implication of ([8]) and (Q is that whenever 3^(2;) < there is a solution ipk 
with k = iz/2 that belongs to L^(]R), i.e., the spectrum consists of the obvious real, nonnegative, 
continuous part and a single eigenvalue: E := —z'^/i. 
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The presence of a pair (v4_, of arbitrary constants in the expression for the eigenf unctions 
ijjk is an indication that the energy levels are doubly degenerate. The application of the program of 
pseudo-Hermitian quantum mechanics [2] requires the construction of a complete biorthonormal 
system consisting of the eigenfunctions of H and H'''. The first step in this direction is to make a 
convenient choice for basis eigenfunctions within each degeneracy subspace. It is most convenient 
to choose one of these eigenfunctions reflectionless [23] - |2S]- This would allow one to express 
this eigenfunction using a formula that is valid for both x < and x > 0. It is not difficult to 
see that sin(/cx) is such a reflectionless eigenfunction. Moreover, the fact that sin(fcx) is an odd 
function suggests to choose the second basis eigenfunction to be even0 As we will see below, this 
choice simplifies the imposition of the biorthogonality conditions considerably. 

Denoting the eigenfunctions by ipay "with a = 1, 2 being the degeneracy label, we set 

Vii^(x) := ^ sin(A;x), (10) 

(U) 



1 



-L Z 

■02 (x) := cos(fcx) H sin (fcx) sign (x 



2k 

We can construct the following eigenfunctions of by replacing z by z* in the above formulas. 



X 



X 



^ sin(A;x), (12) 
vr 



1 

TT 



Z 

cos(kx] + — sin(A;x) sign(x) 
k 



(13) 



Next, we wish to check the validity of the biorthonormality relation for the system {ipa^^'a}- 
Using the well-known integral representation of the delta function, 6{k) = (27r)^^ e^^^dx, the 
parity of the eigenfunctions, and the fact that A; + g > 0, we immediately find 

= S{k - q), i^P'M) = i^P'M) = 0. (14) 

Much more complicated is the derivation o^ 

imi) = il + ^)5ik-q) + nz*5ik)5iq) 

= {l + ^)6{k-q), for k,q>0. (15) 

As seen from ( 1T5|) . the system {ipa^^'a} f^i^s to be biorthogonal for z = ±2ik, because 02 is 
orthogonal to both ipf and for all q G M"^. This is an indication of the presence of a spectral 
singularity, namely E = —z'^/A, which occurs whenever z is purely imaginarjo and is consequently 
embedded in the continuous spectrum of H. 



"'^This is possible, because the Hamiltonian is parity-invariant ('P-synimetric). 
^Here we used the identity 6{k) ~ Um„^oo sin(nfc)/(7rfc). 
•^This is because k e M.~^ and z — ±2ik. 
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A spectral singularity is a serious defect that rules out the operator as a viable candidate for a 
physical observable. We will therefore only consider non-imaginary couplings z. In particular, we 
will focus our attention on the cases where z (and hence Q has a positive real part {^{z) > 0), so 
that there is no eigenvalue and the spectrum is M"*" U {0}. In this case, we can define 

0f(x) := ^Bin(fa), '=°^(fa) + &^"'(fa)s.g.l(x) ^^^^ 



TT 



TT 



which satisfy the biorthonormality condition 

{^lW^ = b^,b{k-q). (18) 

Clearly, this relation is invariant under the transformations 

^^•(x) - ^:'(x) := iV,(z, fc)^^'(x), 0^(x) - 0l^x) := iV,(z, A;)-^>^(x), (19) 

where A^q : C x M+ ^ C are functions that tend to 1 as — > oo and Na{z, k)^^ do not vanish 
except possibly for imaginary values of z. 



3 Construction of a Metric Operator 

Extending the results of [121 E] to the model under investigation, we wish to construct a positive- 
definite metric operator of the form [5l [26l [27] 

poo 

V+ = J2 (20) 
a=l 

This operator defines a positive-definite inner product (■, ■)+ := {■\r]+-) that renders H self-adjoint 
and specifies the physical Hilbert space TCphys of the model [;2j. Note, however, that the metric 
operator ( |20l) is not unique [28] - [31]; one can use (j)'^ of (|T9l) to construct other admissible metric 
operators. Indeed, the determination of the coefficient functions Na{z, k) that would reproduce the 
usual metric operator (//+ = 1) and the {L'^-) inner product in the Hermitian limit {^{z) ^ 0) is a 
very difficult problem. The only guiding principle is to make a simple choice for the biorthonormal 
system that shares the symmetries of the Hamiltonian. In the following we will see that the choice 
(IT6l)-( |T71) made in the preceding section does indeed fulfill this highly nontrivial requirement. 
Having made a choice for 0^, we can try to compute the integral kernel for 77+, namely 

r/+(x,y) := (x|7;+|y) = dk 0^(x)0^(y)*. (21) 

a=l 
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Substituting f|T6l) and f|T71) in this equation and simplifying the result, we obtain 
^?+(x,y) = ^[(5(x-y) -5(x + y)] +a(x,y) + z/3(x,y) sign(y) + 



where 



^*/?(y,x) sign(x) + |2;p7(x,y) sign(x) sign(y), (22) 



1 cosffcx) cos(/cy) 1 / N ^ / M / N 

a(x,y) := — / dk = [/^(x + y) + /^(x - y)], (23) 

^/ N 1 ,, cos(A;x) sin(A;y) 1 / x :r / /^.x 



/„(r) := / c/A; . ^ = / dk (26) 

r G M, and n = 0, 1,2. As seen from (122!) - (126|1 . the calculation of rj+{x,y) reduces to that of 
In{f^)- The latter cannot be evaluated in a closed form. We will construct a series expansion for 
/„(r) that would allow for a perturbative treatment of the problem. 
First, we introduce 



r:f .:=i4 = ,A (27) 



4 ' 4 ' 3?(2) V a 

^{z)r 2k k 



2 ^ ' ^ ■ gfj(z) v^' 

where 3? and stand for the real and imaginary part of their argument, respectively^ 
It is not difficult to show that 



:28i 



/oo 
dq 
■oo 



(29) 



where 

[e2 + 2(l-g2)]^2N-i/2 



/(g>^)^=^l+ ^ (g^ + 1)^ J • 
In view of these equations, we may use e as an appropriate non-Hermiticity parameter. In the 
following we will construct a perturbative expansion of the metric operator in terms of e. 
Expanding /(g, e) in powers of e, we have 



*Note that 5ft(2;) and a are both positive. 



6 



where O(e^) stands for terms of order and higher in powers of e. Substituting (IHTi) in (12^ and 
evaluating the resulting integrals yields, after a very lengthy calculation that is partly done by 
Mathematica, the following remarkably simple result. 



h{r) 
hir) 



27rv^ <^ 5{s) 



+ 



-1^1 



is' -?>\s\ + 11 



it: sign(sj e 









I- 


^(1-1^1)1^1 









-1^1 



1 + 



[s' + lsl + r 



s=y/a r 
s=y/a r 



(32) 
(33) 
(34) 



Having obtained In{r), we are in a position to derive an explicit perturbative expansion for 
the metric operator: 



Af-l 



(x,y) = 5^r/f)(x,y) 6™ + O(e^) 



(35) 



m=0 



where = 1, 2, 3, ■ ■ ■ and ri^"^\x, y) is independent of e. Inserting 



in ([23]) - ([25]), using 

the resulting expression to write (122]) in the form (l35l) . and employing various properties of "sign", 
particularly 

sign(x + y)[sign(x) + sign(y)] = 1 + sign (x) sign (y) = 2 6'(xy), 

where 6{x) := [1 + sign(x)]/2 is the step function, we find after miraculous cancellations of a large 
number of terms 



v+\^^ y) 
v+\^, y) 

v+\^, y) 



- y) 



-5R{^)|x-y|/2 



+ ^(-xy) e 



-5R(2)|x+y|/2 



] sign(y 



2 x^l 



^ {[-K(z)|x - y| ^^(xy) + ^^(-xy)] e"^^^)!^-^!/^ 



v+\^, y) 



32 



+ 



[-K(^)|x + y| ^^(-xy) + ^^(xy)] e-«(^)l^+^l/2} , 

{e(xy)|x - y|(l - l^{z)\x - y|) e-«(^)l^-^l/^+ 
^(-xy)|x + y|(l - \^{z)\x + y|) e"«(^)l^+^l/2} sign(y 



2 x^l 



(36) 
(37) 

(38) 
(39) 



We should emphasize that according to ( l36i) the metric operator 77+ tends to the identity operator 
in the Hermitian limit: e — > 0. This is by no means a trivial expectation. It is a consequence of 
our choice for the biorthonormal system. Moreover rj^^^ (x, y) satisfy the Hermiticity condition, 
?7^^(x, y)* = r7^''(y, x), manifestly. A more important property of ri^^\x,Y) is that they define 
bounded (integral) operators r/^^ in all of L^(]R). This can be established using the fact that the 
integrals |?7^^(x, y)|(iy are bounded for all x G M, [32l §111.2.1]. Alternatively, we may employ 
the following direct proof of the boundedness of rj^'' for m = 0,1,2,3, |33]. First we observe 
that |?7^''(x, y)| viewed as a function of x has an upper bound, /i^'"^(y), depending on y such that 
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I^oo f^^"^Ky)d'J < oo. This implies that for all E 



rj^^ip IP = / dx 



oo 



dy r]^^\x, y)?A(y) 



oo 



2 



"OO / /"OO 

< / dx 



dj y/^M(x)^M(y) |^(y)| 

poo poo 

< I rfx^(™)(x) / rfy/x(™)(y) / rfx|V^(y)r = c2||V^f, (40) 



— oo \J— oo 

oo POO 



where we have used the identities: \ri^\x,Y)\ < \r]^_^\x,j)\ = \ri^\j,x)\ < /i*^™')(x), and 

the Schwarz inequahty [3l]. This completes the proof of the boundedness of r^^'^ for m = 0, 1, 2, 3. 
It implies that at least up to the third order terms in e, 77+ is a bounded operator acting in all of 

We conclude this section by noting that although the Hamiltonian H is manifestly non-PT- 
symmetric, we can follow the approach pursued in [27] to define an antilinear symmetry generator 
(generalized PT-operator [29]) for this Hamiltonian according to 

a=l 

where ^ is complex-conjugation operator defined by (*('C|)|C) •= (CIO5 '4'iC ^ L^(M). In 

view of (HE]), (HZD, and dS]), 

2 poo poo 

(x|r|y) = Y. dk (-l)>.^x)0^(y) = J2 dk (-l)>,^x)V^„^(y)*. 

a=l „=1 Jo 

We will not attempt to obtain a more explicit expression for r, because unlike the metric operator 
r does not enter in the calculation of the physically relevant quantities. 



4 Equivalent Hermitian Hamiltonian 

Having obtained a perturbative expansion for the metric operator we may proceed with the cal- 
culation of the equivalent Hermitian Hamiltonian [9l [21 HJ [5|, [35] - [39j: 

h := 7]'J^Hr]-'^\ (42) 

Using the exponential representation of the metric operator [ID]: rj^ = with Q = Ylm=i Qm^"^, 
the pseudo-Hermiticity relation [12]: W = rj^Hrj^^, and the fact that Qi = —'r]^^\ we first 
calculate the equivalent Hermitian Hamiltonian corresponding to the dimensionless Hamiltonian 
H. The result is [30] 

h = --^ + ^{z) Six) + h(2) e2 + oie'), (43) 
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where 



h'^':=^hi-".«(x)l. (44) 



In view of f l37jl and (jH]), we can easily compute 

(x|h(^)|y) = ^[<5(y)-5(x)]ry«(x,y), 

= ^ e-«Wlyl/2 + 5(y) e-«(^)l-l/2] , (45) 

where the latter expression is to be treated in the sense of distributions (it is valid inside an 
integral). 

Next, we express the Hermit ian Hamiltonian h in terms of the original (unsealed) physical 
variables and the relevant length scale of the problem which, as we explain below, is given by 



This yields 



h = ^ + 3?(C) Six) + ^iCr h^'^ + OmCf), (47) 
2m 



where h^'^^ is defined in terms of its integral kernel. 



h^'\x, y) := ^ [5{x) e-l^l/^ + 5{y) e'l^'l/^ , (48) 
according to 

poo 

(/i(2V)(x):= / h^^\x,y)ij{y)dy = c^[^]e-\''\"' + ci[ij]5{x), (49) 



,1,1 

Clearly h^"^^ and consequently h are nonlocal operators |2]. Furthermore, they are both real (T- 
symmetric) and P-symmetric|f] 

To demonstrate the physical consequences of the imaginary part of the coupling constant C, 
in the original Hamiltonian ([T]) and appreciate the meaning of the length scale L, we calculate 
the energy expectation value for a Gaussian position wave function centered at x = and having 
mean momentum p = hk and width a, 

^(x) = (7ra2)-i/4 e-^+''=". (50) 

It is important to note that we work solely in the Hermitian representation of the quantum system 
where h represents the Hamiltonian of the system and x the position operator. Naturally, we view 
as an element of Ti which yields the probability density of the localization of the particle in the 
physical space as |\E'(x)p. The corresponding element of T^phys is given hy ip = r]^^^'^'^ , [2]. 



^Another property of /i*^^^ is that it eliminates odd wave functions. 
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Figure 1: Plots of as a function of cr for A; = (the full curve), k = 1 (the 
dotted curve), k = 2 (the dashed curve), and k = 4 (the dashed-dotted curve) 
in units where L = 1. 

The energy expectation value of a particle in the state described by the (normalized) position 
wave function flSOj) has the form 



where 



Q{a, k) := e 



h\a-'^ + 2k'^) 3fJ(C) fmn{a,k) 



Am 



+ 



23/2^2 



■i{fc2-L-2)a2 



cos{L~^ka') - 3? \ e*^"^'^"^ erf [2^1/2(^-1 ^ 



(51) 



(52) 



and erf (x) := 27r^^/^ J^^ e~^^dy is the error function. The presence of the exponential factor on the 
right-hand side of fl52|) suggests that the non-Hermiticity effect decays rapidly for mean momentum 
values p = hk outside the range [—hL~^, hL~^]. Figure 1 shows the plots of f2 as a function of a 
for various values of k. As one increases \k\ the non-Hermiticity effect diminishes. The maximum 
is attained for a stationary Gaussian wave packet {k = 0) for which 



,,,,,,, 5J(C) / me^ ri-erf(2-i/2L-i 



^icr+omcY). (53) 



23/2 /j2 



It is not difficult to show that for such a wave packet 

^ ^ + 5S + {j^) 3(0^ + o ((§)3) + omf) 



— ^2- + + 



1 



2 /'Z'i _l_ 



for a ^ L, 
9(C)2 + 0((f)3) +0(9(C)3) for cT<CL. 



Next, we compute the energy expectation value for a stationary Gaussian wave packet of width 
a and mean position a, 

^(x) = (Tra^)-!/^ e"^. (54) 
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Figure 2: Plots of F as a function of a for a = 0.5 (the full curve), a = 1 (the 
dotted curve), cr = 2 (the dashed curve), and a = 3 (the dashed-dotted curve) 
in units where L = 1. 



The result is 



Araa'' 



+ 



_ «2 

e 



TT a 



mF(cr, a] 
23/2/^2 



(55) 



where 



Ffo", a) 



2 2 \ 

1 a 1 

7^ J 



{cosh(f ) - leierf[2-V2(| + ^)] _ le'f erf [2-V2(z _ ^)]} 



Figure 2 shows the plots of F as a function of a for various values of a. As seen from these plots 
the non-Hermiticity effect is substantially smaller for mean positions outside [—L,L]. Therefore, 
L determines the range of the non-Hermitian (nonlocal) interaction. 

In order to make a crude estimate for the magnitude L, consider the application of the delta 
function potential in modelling a point defect in a one-dimensional electron gas system. If we take 
the spatial size d of the defect (lattice size of the crystal) to be of the order of 1 Angstrom and 
the strength of the real part of the potential^ to be of the order of 1 ev, for an electron (of usual 
mass) we find L to be of the order 10^^° Angstrom! Similarly we can obtain an order of magnitude 
estimate for the strength of the non-Hermitian interaction namely m^lz)"^/ (Sh'^d). This turns out 
to be lO^e^ ev. Therefore, to ensure the validity of our perturbative calculation of h, we need to 
take e <^ 10~^. We also recall that the non-Hermitian interaction will be significant, if it is stronger 
than the thermal effects, i.e., m^{zY / {ifv^d) > kT. At room temperature {kT ^ 10^^ ev), this 
implies e > 10~^. 

Finally, we would like to point out that the calculation of the energy expectation value can 
be performed using the pseudo-Hermitian representation of the system. This requires calculation 

— 1/2 

of the state vector ip = ri_^_ corresponding to the position wave function \E'. The energy 



^Here we write the real part of the potential as d ^di{z)d{d ^x) and identify d ^di{z) with its strength. 
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expectation value then takes the form {tpjHijj)^ = {ijj\ri+Hip). This calculation is by no means 
easier to perform than the one reported above [Ullll]. It simplifies to some extend, if one chooses 
ijj directly, e.g., identify iplx) with a Gaussian wave packet. However, note that iplx) is void of a 
direct physical meaning; it is not the position wave function for the state it describes. In order 
to assign a physical meaning for ip in terms of the position of the particle, one must compute the 

1/2 

corresponding position wave function [2j, namely = i]^ ip. 

5 Conclusion 

In this paper we applied the machinery of pseudo-Hermitian quantum mechanics to explore a 
unitary quantum system determined by a delta function potential with a complex coupling con- 
stant (. For an imaginary coupling constant there exists a spectral singularity. For 9ft(C) > 
the spectrum is purely continuous and one can construct a complete biorthonormal system. The 
double degeneracy of the spectrum complicates the choice of a biorthonormal system. We selected 
an appropriate biorthonormal system that simplified the calculations and had a symmetric ex- 
pression for the pair of basis eigenfunctions associated with each degeneracy subspace. We then 
constructed the corresponding metric operator 77+ perturbatively and showed that it tended to 
the identity operator in the non-Hermitian limit and was indeed a bounded operator at least up 
to and including the third order terms that we computed. This is quite remarkable, for there are 
an infinity of other biorthonormal systems such that the corresponding metric operator is either 
unbounded or fails to yield the identity operator in the Hermitian limit. 

Next, we constructed the equivalent nonlocal Hermitian Hamiltonian h for the system. The na- 
ture of the nonlocality of h is quite intriguing, because it originates from a complex delta-function 
potential which is actually ultra-local! This seems to be the reason why the non- Hermit icity effect 
appears in the Hermitian Hamiltonian in the form of a short range interaction, i.e. it decays 
rapidly outside the interaction region: [— L,L], where L = ^7^^- To establish this we calculated 
the expectation value of energy for various Gaussian position wave functions. For a non-stationary 
Gaussian wave packet centered at the origin, the non-Hermitian effect reaches its maximum for 
mean momenta in the range [—hL^^,hL^^]. For a stationary wave packet, it becomes sizable 
whenever the mean position of the packet lies within the interaction region [—L, L]. 

The results reported above show how the methods developed for treating systems with a 
discrete spectrum [3 [2] generalize to specific models with a continuous spectrum. Such a gen- 
eralization has previously been employed in the treatment of the PT-symmetric potential ([3]) 
as reported in [5]. The delta function potential considered in the present paper is manifestly 
non-PT-symmetric, yet we could successfully apply the methods of pseudo-Hermitian quantum 
mechanics [21 El HI [30] to reveal its physical content. 
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